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Abstract. Empirical studies of the first generation of stars and quasars in the Universe 
will likely become feasible over the next decade. The Next Generation Space Telescope 
will provide direct imaging and photometry of sub-galactic objects at 2; ^ 10, while 
microwave anisotropy experiments, such as MAP or Planck, will set constraints on the 
ionization history of the intergalactic medium due to these sources. We describe the 
expected signals that will be detectable with these future instruments. 



INTRODUCTION 

The cosmic dark age [1] ended when the first gas clouds condensed out of the 
primordial fluctuations at redshifts z = 10 — 20 [2] . These condensations gave birth 
to the first star-clusters or mini-quasars in the Universe. Previous observations 
have not yet probed this epoch. Existing optical or infrared telescopes are capable of 
reaching out to redshifts z ;^ 6, while current anisotropy experiments of the cosmic 
microwave background (CMB) probe only the recombination epoch at 2 ~ 10^. The 
remaining observational gap between these redshift regimes is likely to be bridged 
over the coming decade. In particular, the planned launch of the infrared Next 
Generation Space Telescope (NGST) will enable direct imaging of sub-galactic 
objects at z ^ 10, while microwave anisotropy satellites such as MAP or Planck 
will measure complementary signatures from the reionization of the intergalactic 
medium at these redshifts. 

The detection [3,4] of Lya emission from galaxies at redshifts up to 2; = 5.7 
demonstrates that reionization due to the first generation of sources must have 
occurred at yet higher redshifts; otherwise, the damping wing of Lya absorption 
by the neutral intergalactic medium would have eliminated the Lya line in the 
observed spectrum of these sources [5,6]. In this review, we predict the redshift 
of reionization in popular cosmological models, and estimate the signals from this 
epoch that will become detectable with the NGST, MAP, or the Planck satellites. 
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PROPERTIES OF THE FIRST OBJECTS 

The first objects left beliind a variety of fossil evidence of their existence [7], in- 
cluding the enrichment of the intergalactic medium (IGM) with heavy elements, the 
reionization of the IGM, the distortion of the CMB spectrum by dust-processed ra- 
diation [8], and the production of stellar remnants. In the following, we explore the 
signatures of the early stars and mini-quasars using a simple semi-analytical model, 
based on the Press-Schechter formalism [9]. We calibrate the total amount of light 
that stars or mini-quasars produce based on data from redshifts z ^ 5. The effi- 
ciency of early star formation is calibrated based on the observed metallicity of the 
intergalactic medium [10,11], while the early quasars are constrained so as to match 
the quasar luminosity function [12] at redshifts 2; ;^ 5, as well as data from the Hub- 
ble Deep Field (HDF) on faint point-sources. We focus on a particular cosmological 
model with a cosmological constant, (Qq, ^a, ^h, h, Ug/i-i, n)=(0.35, 0.65, 0.04, 0.65, 
0.87, 0.96), named the "concordance model" by Ostriker & Steinhardt [13]. For a 
discussion of other cosmological models, as well as a more detailed description of 
our methods and results, we refer the reader to several papers [26,14-20]. More ad- 
vanced 3-D numerical simulations have only now started to address the complicated 
physics associated with the fragmentation, chemistry, and radiative transfer of the 
primordial molecular clouds [21], as well as with the reionization of the IGM [22]. 

Following collapse, the gas in the first baryonic condensations is virialized by a 
strong shock [23] . The shock-heated gas can only continue to collapse and fragment 
if it cools on a timescale shorter than the Hubble time. In the metal-poor primordial 
gas, the only coolants that satisfy this requirement [24] are neutral atomic hydrogen 
(H) and molecular hydrogen (H2). However, H2 molecules are fragile, and are easily 
photo-dissociated throughout the universe by trace amounts of starlight [25,26] 
that are well below the level required for complete reionization of the Universe. 
Hence, most of the sources that ionized the Universe formed inside objects with 
virial temperatures Tvu ^ lO^K, or masses ~ 10^ M©, which cooled via atomic 
transitions. Depending on the details of their cooling and angular momentum 
transport, the gas in these objects fragmented into stars, or formed a central black 
hole exhibiting quasar activity. Although the first objects contained only a small 
fraction of the total mass of the universe, they could have had a dramatic effect 
on the subsequent evolution of the ionization and temperature of the rest of the 
gas [27]. Since nuclear fusion releases ~ 7 MeV per baryon, and accretion onto a 
black hole may release even more energy, and since the ionization of a hydrogen 
atom requires only 13.6 eV, it is sufficient to convert a small fraction of the baryonic 
mass into either stars or black holes in order to ionize the rest of the Universe. 

The cooling gas clouds eventually fragment into stars [28]. Although the actual 
fragmentation process is likely to be rather complex, the average fraction /star of the 
collapsed gas converted into stars can be calibrated empirically so as to reproduce 
the average metaUicity observed in the Universe at 2; f« 3. For the purpose of this 
calibration, we use the average C/H ratio, inferred from CIV absorption lines in 
Lya forest clouds [10,11]. The observed ratio is between 10~^ and 10~^ of the solar 
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value [29] . If the carbon produced in the early star clusters is uniformly mixed with 
the rest of the baryons in the Universe, this implies /star ~2-20% for a Scalo [30] 

stellar mass function. This number assumes inefficient hot bottom burning, i.e. 
maximal carbon yields [31], and includes a factor of ~3 due to the finite time 
required to produce carbon inside the stars (in a Press-Schechter star formation 
history, only a third of the total stellar carbon yield is produced and ejected by 
2; = 3). Ultimately, 3-D simulations of the first generation of stars might be used 
to infer the expected star-formation efficiency in the first generation of gas clouds. 
Preliminary runs [32] imply that ~ 1% of the gas condenses into dense cores which 
could yield massive stars. 

An even smaller fraction of the cooling gas might condense at the center of the 
potential well of each cloud and form a massive black hole, exhibiting mini-quasar 
activity. In the simplest scenario we postulate that the peak luminosity of each 
black hole is proportional to its mass, and there exists a universal quasar light- 
curve in Eddington units. This hypothesis is motivated by the fact that for a 
sufficiently high fueling rate, quasars are likely to shine at their maximum possible 
luminosity, which is some constant fraction of the Eddington limit, for a time which 
is dictated by their final mass and radiative efficiency. Allowing the final black hole 
mass Mbh to be a fixed fraction of the total halo mass Mhaio, we find that there 
exists a universal light curve [L{t) — LEddexp(— t/to), with to ~ 10^ yr], for which 
the Press-Schechter theory provides an excellent fit to the observed evolution of the 
luminosity function (LF) of bright quasars between redshifts 2.6 < z < 4.5. The 
required black hole to halo gas mass ratio is Mbh/Mgas ~ lO^^-^Qo/^b = 5.5 x 10~^, 
close to the typical value of ~ 6 x 10~^ found for the ratio of black hole mass to 
spheroid mass in a dozen nearby galaxies [33,34]. The existence of massive black 
holes in the centers of low-mass galaxies such as M32 or NGC 4486B [33,34] implies 
that the process of black hole formation does not discriminate against galaxies of 
this type. Since galaxies at 2; ~ 10 have a similar mass and velocity dispersion 
as these low-redshift examples, it is conceivable that low-luminosity quasars con- 
tributed significantly to the reionization of the Universe. 

One does expect, however, that the ratio Mbh/Mgas would have a substantial in- 
trinsic scatter. Observationally, the scatter around the average value of \og{My^^/ L) 
is 0.3 [34], while the standard deviation in log Mbh/Mgas has been found to be 
0" ~ 0.5 [35]. Such an intrinsic scatter would fiatten the predicted quasar LF at 
the high mass end, where the LF is a steep function of black hole mass. As an 
illustrative example, we show in Figure 1 the mass function of black holes in the 
Press-Schechter model of halos, with or without this scatter (solid versus short- 
dashed lines). In order to eliminate the fiattening introduced by the scatter, we 
find that the average black hole to halo mass ratio must be reduced by ~ 50%. 
The dot-dashed line in Figure 1 demonstrates that such a reduction would in- 
deed compensate for the effect of an intrinsic scatter in the relevant mass range 
(10^ Mq ^ Mbh ^ 10^° Mq). Figure 1 also shows the effect of a more significant 
intrinsic scatter (cr ~ 1) on the black hole mass function (long-dashed lines). We 
find that the predicted black hole mass function in the presence of such a large 
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FIGURE 1. The comoving mass function of black holes at redshift z = 2> when a scatter is 
introduced to the logarithm of the ratio between the black hole and halo masses, log(Mbh/Afhaio)- 

scatter would be significantly different from any model with a constant value for 

In reality, the relation between the black hole and halo masses may be more 
complicated than linear in reality. With the introduction of additional free param- 
eters, a non-linear (mass and redshift dependent) relation between the black-hole 
and halo masses can also lead to acceptable fits [36] of the observed quasar LF. 
The nonlinearity in the relation must be related to the physics of the formation 
process of low-luminosity quasars, which was discussed in several papers [36-38]. 
If the black hole formation efficiency decreases in smaller halos, this would flatten 
the faint end of the LF, and therefore could not compensate the effect of a large 
intrinsic scatter of the type shown in Figure 1. Indeed, in order to fit the bright 
end of the LF in a model with a large intrinsic scatter, one must postulate that the 
black hole formation efficiency decreases in larger halos. 



INFRARED NUMBER COUNTS 

The Next Generation Space Telescope {NGST, [39]) will be able to detect the 
early population of star clusters and mini-quasars. NGST is scheduled for launch 
in 2007, and is expected to reach an imaging sensitivity of ~ 1 nJy (S/N=10 at 
spectral resolution A/ A A = 3) for extended sources after several hours of integration 
in the wavelength range of l-3.5;um. Figure 2 shows the predicted number counts in 
the models described above, normalized to a 5' x 5' field of view. This figure shows 
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FIGURE 2. Infrared Number Counts. The solid curves refer to quasars, while the long/short 
dashed curves correspond to star clusters with low/high normalization for the star formation 
efficiency. The curves labeled "5" or "10" show the cumulative number of objects with redshifts 
above z = 5 or 10. 

separately the number per logarithmic flux interval of all objects with redshifts 
z > 5 (thin lines), and 2; > 10 (thick lines). The number of detectable sources 
is high; NGST will be able to probe about ~ 100 quasars at 2; > 10, and ~ 200 
quasars at ^ > 5 per field of view. The bright-end tail of the number counts 
approximately follows a power law, with dN/dFy oc F~'^'^. The dashed lines show 
the corresponding number counts of "star-clusters" , assuming that each halo shines 
due to a starburst that converts a fraction of 2% (long-dashed) or 20% (short- 
dashed) of the gas into stars. These fines indicate that NGST would detect ~ 
40 — 300 star-clusters at z > 10 per field of view, and ~ 600 — 10^ clusters at 
z > '0. Unlike quasars, star clusters could in principle be resolved if they extend 
over a scale comparable to the virial radius of their dark matter halos [16]. The 
supernovae and 7-ray bursts in these star clusters might outshine their hosts and 
may also be directly observable [40,41]. 



CONSTRAINTS FROM THE HUBBLE DEEP FIELD 



High resolution, deep imaging surveys can be used to set important constraints 
on semi-analytical models of the type described above. The properties of faint 
extended sources found in the Hubble Deep Field (HDF) [42] agree with detailed 
semi-analytic models of galaxy formation [43]. On the other hand, the HDF has 
revealed only a handful of faint unresolved sources, but none with the colors ex- 
pected for high redshift quasars [44]. The simplest mini-quasar model described 
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FIGURE 3. V, I, J, and H counts for mini-quasars in a model with a minimum halo circular 
velocity of Vdic = 75 km s~^. This model is consistent with HDF data in the V band at the 5% 
level. Sensitivities are shown in each hand for the same signal-to-noise ratio and exposure time 
as the optical HDF. 

above predicts the existence of ~ 10 B-band "dropouts" in the HDF, inconsis- 
tently with the lack of detection of such dropouts up to the ~ 50% completeness 
limit at ~ 29 in the HDF. To reconcile the models with the data, a mechanism 
is needed for suppressing the formation of quasars in halos with circular velocities 
Vcuc ^50 — 75 km s~^. This suppression naturally arises due to the photo-ionization 
heating of the intergalactic gas by the UV background after reionization [45,46]. 
Alternative effects could help reduce the quasar number counts, such as a change in 
the background cosmology, a shift in the "big blue bump" component of the quasar 
spectrum to higher energies due to the lower black hole masses in mini-quasars, or 
a nonlinear black hole to halo mass relation; however, these effects are too small to 
account for the lack of detections in the HDF [20] . 

The mini-quasars may not necessarily appear as point sources in the HDF if 
their extended host galaxies are actually resolved by HST. In fact, twelve candidate 
sources of activity in the nuclei of galaxies at high-redshifts {z > 3.5) have recently 
been identified [47] in the HDF. All of these point-like sources are embedded in 
extended host galaxies which are relatively bright (l^ ~ 26 — 27) and outshine their 
AGNs by typically 1 mag. As a consequence, these AGNs would have been missed 
by previous searches for isolated point-like sources [44]. In the models described 
above, the mini-quasars peak at a flux ~ 1 — 3.5 mag brighter in V than their 
host galaxy, which is assumed to undergo a starburst inside the same halo [18]. 
The twelve HDF candidates must reflect faint AGN activity in bright galaxies of 
relatively massive halos (analogous to a weak Seyfert activity), rather than faint 
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AGN activity in small halos as expected in our model. The observed sources might 
imply a phase in the history of massive halos that corresponds to an additional low- 
luminosity tail of the quasar lightcurve discussed above. An extended lightcurve 
would still be consistent with the luminosity function derived from the bright AGN 
phase and could explain the existence of the embedded AGNs in the HDF. The 
detection of faint embedded AGNs could also be related to the intrinsic scatter 
in the distribution of Mbh/Mhaio relation, and reflect objects with unusually small 
values of this mass ratio. 

The longer infrared wavelengths, such as the J and H infrared bands, are better 
suited for studying the Universe at z ^ 5. Forthcoming data on point-sources from 
NICMOS observations of the HDF [48] could improve the constraints on mini- 
quasar models. In Figure 3, we show the expected number counts of mini-quasars 
in the V , I. J. and H bands for the model which is consistent with the optical 
HDF data. The number of objects predicted in the /, J, and H bands is higher 
than in the optical HDF. The NICMOS data in these bands would either reveal 
several high redshift mini-quasars or else place tighter constraints on quasar models 
than currently possible using V and / data. With the post-reionization feedback 
on halos imposed by the optical HDF data, Vdrc ^ 75 km s~^, we still expect at 
least ~ 5 mini-quasars to be found at z > 3.5 in the NICMOS J and H bands. 
A non-detection by NICMOS would translate to a minimum circular velocity of 
^circ ^100 km s~^, or a factor of ~ 2 increase in the low-mass cutoff for halos 
harboring quasars. 

WAS THE UNIVERSE REIONIZED BY STARS OR 

MINI-QUASARS? 

Given either the star-formation or quasar black-hole formation histories, we de- 
rive the reionization history of the IGM by following the radius of the expanding 
Stromgren sphere around each source. The reionization history depends on the 
time-dependent production rate of ionizing photons, their escape fraction, and the 
recombination rate of the IGM, which are all functions of redshift. The production 
rate of ionizing photons per quasar follows from the median quasar spectrum [49] 
and the light-curve we derived. The analogous rate per star follows from the time- 
dependent composite stellar spectrum, constructed from standard stellar atmo- 
sphere atlases [50] and evolutionary tracks [51]. We computed the escape fraction 
of ionizing photons in each halo, assuming ionization equilibrium inside an isother- 
mal sphere where both the stars and the gas are distributed with a profile. For 
quasars, we assumed that the escape fraction is 100%. During reionization, we also 
assumed that the formation of low-mass halos is suppressed by photo-ionization 
heating inside the already ionized cosmological HII regions. 

Figure 4 summarizes the resulting reionization histories from stars or quasars 
in our models with a ACDM cosmology. The results for stars are shown in two 
cases, one with Ziqm — 10~^Zq (dashed lines) and the other with Ziqm — 1O~^Z0 



FIRST STARS AND QUASARS 




FIGURE 4. Reionization history. Clockwise, the different panels show: (i) the collapsed fraction 
of baryons; (ii) the background flux at the Lyman limit; (Hi) the volume filling factor of ionized 
hydrogen; and (iv) the optical depth to electron scattering, and the corresponding damping factor 

for the power spectrum decomposition of microwave anisotropics as a function of the spherical 
harmonic index I. The solid curves refer to quasars, while the dotted or dashed curves correspond 
to stars with a low (2%) or high (20%) normalization for the star formation efficiency. 



(dotted lines), to bracket the allowed IGM metallicity range. The panels in Figure 4 
show (clockwise) the total collapsed fraction of baryons available for star or quasar 
formation; the evolution of the average comoving flux, J21 at the local Lyman 
limit frequency, in units of 10~^^ erg s~^ cm~^ Hz^-*- sr^"^; the resulting cvohition 
of the ionized fraction of hydrogen, -Fhh; and the consequent damping of the CMB 
anisotropies. In the first two panels, we imposed Vcuc > 75 km s~^ at z < 2;reion- 
The dashed and dotted curves indicate that stars ionize the IGM by a redshift 
9 ^ z ^ 13; while the solid curve shows that quasars reionize the IGM at z ^ 11. 
This result can be understood in terms of the total number of ionizing photons 
produced per unit halo mass; given our normalizations of the efficiencies of star 
and quasar black hole formation, the relative ratios of this number in the three 
cases are 1 0.37 0.1, respectively. A comparison of our quasar and stellar 
template spectra shows that stars will not reionize Hell, while quasars reionize 
Hell at essentially the H reionization redshift. Therefore, recent claims that HcII 
reionization might have been observed at 2; ~ 3 [52] could rule-out the presence 
of mini-quasars with hard spectra extending to X-rays at high redshifts, if these 
claims are verified by future observations. 

The X-ray background (XRB) [53,54] might provide another useful constraint 

on the mini-quasar models. In Figure 5 we show the predicted spectrum of the UV 
to the soft XRB at ^ = in these models (solid lines). In computing the spectrum. 
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FIGURE 5. Spectrum of the UV/X-ray background in the mini-quasar model, assuming that 
the median X-ray spectrum of quasars [49j is universal. The solid curves show the spectrum with 
and without absorption by the Lya forest at high-redshifts. The short and long dashed lines show 
the unresolved fraction (assumed to be 25%) of the observed X-ray background spectrum (from 
[53] and [54j). 

we included absorption by neutral H and He in the IGM at ^ > ^reion = 11 and 
hydrogen absorption [55] by the (extrapolated) Lya forest at z < 11. Also shown 
in this figure is the unresolved 25% fraction of the observed soft XRB [53,54]. The 
dashed lines in figure 5 represent an upper limit on any component of the XRB 
that could arise from high-rcdshift quasar activity. As the figure shows, the mini- 
quasar models overpredict the unresolved flux by a factor of ~ 2 — 7 in the 0.1-1 keV 
range, as they produce a flux comparable to the entire soft XRB flux. If an even 
larger fraction of the XRB will be resolved into low-rcdshift AGNs in the future, 
then the XRB could be used to place stringent constraints on the X-ray spectrum 
or the abundance of the mini-quasars discussed here. 

CAN THE REIONIZATION REDSHIFT BE INFERRED 
FROM A SOURCE SPECTRUM? 

The spectrum of a source at a rcdshift Zs > -Zrcion should show a Gunn-Peterson 
(GP) [56] trough due to absorption by the neutral IGM at wavelengths shorter than 
the local Lya resonance at the source, Aobs < ^a{^ + Zg)- By itself, the detection of 
such a trough would not uniquely establish the fact that the source is located beyond 
-^reion, siuce the lack of any observed flux could be equally caused by: (i) ionized 
regions with some residual neutral fraction, (ii) individual damped Lya absorbers, 
or (iii) line blanketing from lower column density Lya forest absorbers. On the 
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FIGURE 6. Spectrum of a source at Zg = 7.08, assuming sudden reionization at a redshift 
-Jreion = 7. Thc soUd cuTvcs show the spectrum without absorption by the high-redshift Lya forest, 
and the dashed lines show the spectrum when the damping wings are also ignored. 



other hand, for a source located at a redshift Zg beyond but close to reionization, 
(1 + -Zrcion) < (1 + -^s) < |f (1 + -Zrcion), the GP trough splits iuto disjoint Lyman 
a, P, and possibly higher Lyman series troughs, with some transmitted flux in 
between these troughs. Although the transmitted flux is suppressed considerably 
by the dense Lya forest after reionization, it is still detectable for sufficiently bright 
sources, and can be used to infer the reionization redshift. 

As an example, we show in Figure 6 the simulated spectrum around the Lyman 
a and f3 GP troughs of a source at redshift Zs = 7.08, assuming that reionization 
occurs suddenly at 2;reion = 7. We have included the extrapolated effects of Lya 
absorbers along the hnes of sight, whose statistics were chosen so as to obey the 
redshift dependence and absorption line characteristics of observational data at 
z < 4.3 [57]. Although the continuum flux is strongly suppressed, the spectrum 
contains numerous transmission features; these features are typically a few A wide, 
have a central intensity of a few percent of the underlying continuum, and are 
separated by ~ lOA. For sudden reionization, the nominal integration time of 
about 10 hours for a ~ 10 nJy sensitivity (with A/AA=100 and S/N=10) would 
be sufficient to determine ^rcion up to a redshift of ~ 7 with a high precision, 
by determining the location of the short-wavelength edge of the troughs shown in 
Figure 6. More gradual reionization would smear the edge of the GP trough. Based 
on the extrapolation of the Lya forest to higher redshifts, the required sensitivity 
needs to be one or two orders of magnitude higher if ^^reion ~ 8 or ~ 9 (see [19] for 
details) . 
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FIGURE 7. Effect of dust on the background flux. The top panel shows the comoving spectra 
in four different models at z = 3, and the bottom panel shows the corresponding evolution of the 
dust temperature and Compton y-parameter (see text). 



SIGNATURES IMPRINTED ON THE CMB 

The free electrons produced by the reionization of the intergalactic medium par- 
tially smooth-out the temperature anisotropics of the CMB via Thomson scat- 
tering. Given the ionized fraction of hydrogen as a function of redshift, one can 
readily derive the electron scattering optical depth {r^s), as well as the anisotropy 
damping factor (R^), as functions of the spherical harmonic index £ of the multipole 
expansion of the anisotropics on the sky [58] . As illustrated in the lower right panel 
of Figure 4, the amplitude of the anisotropics is reduced by ~ 6 — 10% on scales 
below the angular scale of the horizon at reionization {£ ^ 10). Although small, 
this reduction is within the proposed sensitivities of the future MAP and Planck 
satellites, provided that both temperature and polarization anisotropy data will be 
gathered in these experiments (see Table 2 in [59]). 

In addition, the dust that is inevitably produced by the first type II supcrnovac, 
absorbs the UV emission from early stars and quasars and re-cmits this energy at 
longer wavelengths, where it distorts the CMB spectrum. We have calculated this 
spectral distortion assuming that each type II supernova yields O.SM© of dust with 
the wavelength-dependent opacity of Galactic dust [60]. We have conservatively 
assumed that similarly to the observed intergalactic mixing of metals, this dust 
gets uniformly distributed throughout the intergalactic medium, dumpiness of the 
dust around UV sources would only enhance our predicted spectral distortion. The 
top panel of Figure 7 shows the resulting total comoving spectrum of the radiation 
background (CMB + direct quasar and/or stellar emission + dust emission) at 
z = 3. More distortion could be added between < 2; < 3 by dust and radiation 



FIRST STARS AND QUASARS 



from galaxies. For reference, we also show by the dot-dashed hnes the recently 
detected cosmic infrared background (CIB, [61]), and the typical dust peak in our 
calculations (dotted lines). 

Figure 7 shows that in our models the dust emission peaks at a wavelength which 
is an order of magnitude longer than that of the CIB peak. This is a result of 
our assumption of a homogeneous dust distribution, and the consequent cold dust 
temperature. An inhomogeneous distribution of dust would raise its temperature, 
and could contribute significantly to the observed CIB peak. The deviation from 
the pure 2.728(1 + 2;)K blackbody shape is quantified by the Compton y-parameter, 
whose redshift evolution is shown in the bottom panel. Ignoring the UV flux from 
quasars, we obtain 1.1 x 10~^ < yc < 8.2 x 10~^ at ^ = 3 (dashed hnes), just 
below the upper limit [62] set by COBE, y < 1.5 x 10~^. Adding the UV flux of 
quasars increases the y-paramctcr to 4.1 x 10^^ < yc < 2 x 10^^. The distortion 
by the intergalactic dust may have been overestimated in our calculation as we 
ignored the absorption of the UV background by the neutral component of the 
IGM. Nevertheless, a substantial fraction (~ 10-50%) of the total y-parameter 
results simply from the direct far-infrared emission by early quasars and could be 
present even in the absence of any intergalactic dust. 
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